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Abstract By means of population synthesis we model the Galactic ensemble of helium stars. It is 
assumed that all helium stars are formed in binaries. Under this assumption, single helium stars are produced 
by the mergers of helium remnants of components of close binaries (mainly, by merging helium white dwarfs) 
and by disruption of binaries with helium components in supernovae explosions. The estimate of the total 
birthrate of helium stars in the Galaxy is 0.043 yr -1 , their total number is estimated as 4 x 10 6 . The 
rate of binarity in the total sample is 76%. We construct a subsample of low-mass (Mh c < 2M Q ) helium 
stars limited by observational selection effects: stellar magnitude (Vn e < 16), ratio of stellar magnitudes of 
components in binaries (Vh c < V comp ), lower limit of the semiamplitude of radial velocity that is necessary 
for discovery of binarity [K m in = 30 km / s) . The parameters of this "observable" sample are in satisfactory 
agreement with the parameters of the observed ensemble of sdB stars. In particular, in the selection-limited 
sample binarity rate is 58%. We analyze the relations between orbital periods and masses of helium stars 
and their companions in systems with different combinations of components. We expect that overwhelming 
majority (~ 90%) of unobserved components in binary sdB stars are white dwarfs, predominantly, carbon- 
oxygen ones. 

^ ; 1 Introduction 

> : 

CN| , The aim of the present paper is investigation of nondegenerate helium stars. Among observed objects they 

are represented by two groups - helium subdwarfs and Wolf-Rayet stars (WR). The interest to helium 
stars is enhanced by the fact that they are related to the stars of another groups: AM CVn stars, R CrB 
stars, double-degenerates (possible progenitors of type la supernovae). The most massive helium stars are 
progenitors of type lb and Ic supernovae. The optical component of one of the most remarkable X-ray 
sources - Cyg X-3 is a helium star. Among stars with spectroscopically determined masses the most massive 
is the close binary WR 20a; both components of this system are helium stars with M ~ 80 Mq llRauw et all 
12004 120051) . 

O Helium subdwarfs are subdivided into three groups: sdB stars that have hydrogen-rich atmospheres, 

sdO stars with atmospheres dominated by helium, and an intermediate group of sdOB stars. The masses of 
helium subdwarfs that are usually derived fro m their location in the "effective temperature log T e g - effective 
gravity log g" diagram, are close to 0.5 Mq l|Heherl . Il986t ISaffer et all Il994h . 

£> '• 

. !_| \ These estimates are consi stent with masses det ermined for components of eclipsing systems that range 

^ ■ from 0.48 M Q to 0.54 Mq fsee lDrechsel et all ll200ll) a nd references therein) and with masses obt ained from 

asteroseismological data: 0.49±0.02M Q (PO0014+067 l|Brassard et al!l200lK 0.47 M (Feige 48 foeed et all 
l2004ln . One has to have in mind, however, that the evolutionary tracks of numerous low-mass (M wc j < 
0.3 Mq) helium white dwarfs in the stage of contraction of their hydrogen envelopes, also cross the region of 
the logT c ff - log g diagram populated by helium subdwarfs. Respectively, some of the stars that are classi- 
fied as subdwarfs may happen to be not helium stars but relatively hot helium white dwarfs. For instance, 
the m ass of HD188112 tha t is classified as sdB is 0.24 Mq only, i. e., it is a precursor of a helium white 
dwarf flHeber et all 1200 3a). As another example of uncertainty one may consider HS2333+3927: position 
of this star in the logT c ff - log g diagram "satis fies" both track s of a 0.29 Mq white dwarf and of a 0.47 Mq 
helium star with a thin hydrogen envelope l)Heber et alll20o"4j) . 

The lower limit to the range of masses of WR stars is ~ 7M Q l|Nugis fc La.mersl 120001) . However, 
according to the modern theory of the evolution of close binary stars, mass spectrum of helium stars has to 
be continuous. It is possible that the WR phenomenon is observed for stars with mass > 7 M & only, because 
for mass-luminosity relation typical for the lower mass stars radial pulsations that m ay result in generation 
of sho ck waves that produce dense enough stellar wind are not excited (see, e. g., IFadevev fc Novikoval 

j2nnaft ). 

There are several channels for production of helium stars. Stars with ZAMS masses exceeding ~ 
50 Mq turn into WR stars due to stellar wind mass loss llContil . Il979l) . It is possible that single stars 
and components of wide binaries of moderate mass may almost completely lose their hydrogen envelopes 
immediately before helium flash or during the latter l|Cox & Salneterl . Il96ll iD'Cruz et all Il99fil) . In close 
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binaries, stars with helium cores and thin hydrogen enve lopes form due to the mass loss after Roche lobe 
overflow llKinnenhahn fc WeigerH . ll967llPa,czvriskiLll967ft . Single heliu m stars may form due to the merger 
of he l ium white dwarfs or nondegenerate helium stars in close binaries l|Tutukov fc YungelsonI Il990: Iben, 
IT^ISaiofc Jeffervlliin^. 

The main scenario s for fo r mation of low-mass close binar i es wer e considered hv llhen fc Tutukovl l|l985l ): 
ITutukov fc YungelsonI lll99ft : llhenl l|l990ft : lHan et alJ fewA \20()ti \. In particular, it has been shown that 
the location of helium remnants of components of close binaries and of the products of their merger and 
their evolutio nary track s in the Hertzspmng-R ussell diagram and in the log T c g — log g diagram well fit 
observations (Iben, 1990; Han et al. .. 2002, 2003). Han et al. have considered in detail the relations between 
different formation channels and have shown that, if selection effects are taken into account, it is possible 
to explain the main parameters of observed helium subdwarfs. 

In the present paper, by means of population synthesis, we construct a model of the ensemble of Galactic 
helium stars and analyze the propertie s of its low-mass subsystem (i .e., of helium subd warfs) . We apply 
the p opulation synthesis code used hv ITutukov fc YungelsonI ll2002l l. At difference to lHan et alJ (2002, 
2003), we compare with observations a subset of mo del stars lim i ted b y stellar magnitude, since such a 
limit exists for the main catalog of hot subdwarfs PG l|Green et all Il98fift . We pay main attention to some 
aspects of binarity of hot subdwarfs, in particular, to the distribution of the pairs of subdwarfs with different 
companions over orbital periods. 

Massive helium stars and the products of their evolution are considered in the next paper of the series 
(Paper II). 

In Sect. El we describe the basic assumptions used in the population synthesis code that are immediately 
relevant to the modeling of formation of helium stars. The main results are presented in Sect. El Discussion 
follows in Sect. 0J 



2 Population synthesis for binary stars 

Both high- and low-mass helium stars are formed when one of the components of a close binary system 
overf lows its Roche lobe in the hydrogen shell burning stage l|Kinnenhahn & WeigertL Il967t IPaczvhskiL 
11 967ft (in the co called "case B" of mass exchange) . In the population synthesis code we have the following 
assumptions relevant to the mass exchange in case B. 

Components of close binaries with ZAMS masses M; < 2.8 M© form degenerate helium dwarfs. Com- 
ponents of close binaries that initially had mass higher than 2.8 M (or accumulated that much mass as 
a result of mass exchange) produce nondegenerate helium stars. The mass of the latter is related to Mj 
by the following, appro ximation to the results of evolutionary computations llTvtukov fc YungelsonI . Il973t 
iThen fc Tntnkolrmslh: 

Mhc/Mq = max(O.O66(M i /M ) 1 - 54 , O.O82(M i /M ) 1 - 4 ). (1) 

Equation (QJ gives the minimum mass of helium st ars close to 0.35 M , in a reasonable agreement with 
results of more detailed evolutionary computations l|Ha,n et all l200nft . We neglect the dependence of the 
mass of helium remnant on the instance of the Roche lobe overflow that is not very significant in the case 
B of mass exchange. 

Based on the results of evolutionary computations llTutukov fc YungelsonI . Il973t llben fc Tutukovl 

1985), we approximated the lifetime of stars in the region of Hertzsprung-Russell diagram populated by 
helium stars as 

T ( 7.15 - 3.71og(M Hc /M Q ) for M Ho < 1 M , 

° S HC \7.15-3.71og(MHc/M ) + 2.23[log(M H c/Af )] 1 - 37 forM Ho >lM . 

Stellar wind mass loss in the main-sequence stage is described by initial-final mass relations l|Vanbeverenl 

11991ft : 

_ fMj - 0.00128M2- 5 for M ; < 30 M , 
f ~ |l.77AfP- 763 for M > 30 M . ^' 

It is assumed that in the core-hydrogen burning stage the stars with ZAMS mass in excess of 50 M 
lose by stellar wind so much matter that they become helium stars immediately after completion of main- 
sequence stage and never fill Roche lobes. Thus, only helium stars with mass > 27 M may be products of 
the evolution of single stars or of components of wide binaries. 
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Mass loss by helium stars with Mh c > 5 Mq is described by th e follow ing approximation to the observed 
mass- loss rates by WR-stars (see compilation hy lNugis fc Lament ll2000l )): 



min(1.38 x 10~ 8 (M H e/M Q 



,2.87 



10- 4 )M Q /yr. 



(4) 



The analytical expression in Eq. <@J is suggested by iNelemans fc va.n den Heuvell fcOOlh . The slope of the 
M H e(M He ) relation changes at M He ~ 22.1 M . 

So me signs of the existence of stellar wind are s uspect e d for several hot subdwarfs only llHeber et all 
2003b); applying the model of radiative stellar wind. Fvrml ll2004h estimated the rate of mass loss as M ~ 
10~ n M Q /yr. Taking into account the lifetime of helium stars one may claim that the only effect that such 
a mass loss may have upon evolution of hot subdwarfs is the loss of the remainders of hydrogen envelope 
and, possibly, change of the chemical composition of the subdwarf atmosphere. 

Helium stars with initial mass M\ — 2.8 -j- 5 M© do not expand after helium exhaustion in their cores 
and become carbon-oxygen white dwarfs remaining compact. If M\ = (5 — 11.6) Mq, helium stars finish 
their evolution as carbon-oxygen white dwarfs after expansion in the heli um-shell burning: stage and repeated 
Roche lobe overflow that is accompanied by an insignificant mass loss l|lhen & TutukovL Il98fll) . 

Stars with ZAMS mass Mi = 11.6 ~ 30 Mq produce neutron stars with Mns = 1.4 Mq. Initially more 
massive stars end their lives as black holes. Formation of black holes is accompanied by the loss of 20% of 
the presupernova mass. It is assumed that supernovae explosions are spherically symmetric, i. e., nascent 
neutron stars and black holes do not get kicks (see detailed discussion in Paper II). 

Stars with Mi < 11.6 Mq that overflow Roche lobe in the AGB stage ("case C" of mass exchange) 
produce white dwarfs; the mass of the latter depends on the instance of Roche lobe overflow. More massive 
stars produce neutron stars or black holes, like in case B of mass exchange. 

The relation between initial mass of stars, their mass at the end of main-sequence stage, mass of helium 
stars produced by them, and mass of the final products of evolution is plotted in Fig. 

We have applied t he following assumpti ons on the character of mass-exchange process, based on results 
of numerical modeling JTutukov et allll982h . If both donor and accretor are main-sequence stars, the donor 
has radiative envelope, and mass ratio of components q — M&/M& > 0.8, mass exchange is conservative and 
proceeds in the time scale of the nuclear evolution of the donor. If 0.5 < q < 0.8, the donor loses mass in 
the thermal time scale tthi, but accretion by the main-sequence secondary is limited by the rate Mt%2, that 
corresponds to the thermal time scale of the latter. The mass that may be accreted by the companion is 
limited by AM = M t hitth\- Variation of the separation of components a is calculated under assumption 
that AM is accreted conservatively, but the rest of the donor envelope with mass AM t = Mi — Mf — AM 
is lost from the system; the latter process is described by the equation 
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' M t 1 



(5) 



where J is angular momentum of the binary, M t = Mi + M 2 , 7 = 1.5. The validity of Eq. © with the 
parameter 7 = 1.5 is based on the results of reconstruct ion of the mass-exchange history in observed c lose 
binaries containing white dwarfs and helium subdwarfs JNelemans et all 120001: iNelemans fc Tout! 120051) . 

When the mass ratio of components q < 0.5 or donor has convective envelope extending down to more 
than 50% of stellar radius, mass loss occurs in dynamical time scale and a common envelope is formed. 
Variation of component's separatio n is computed based on the ba lance of the orbital energy of the binary 
and donor envelope binding energy l|Tutukov fc Yungelsonl Il979al) : 



(Mi + to) (M - Mf) 



a; 



Mf TO Mi TO 



Of 



di 



(6) 



In the Eq. © indexes i and f are related to the initial and final mass and separation of components, to 
is accretor mass, a cc is so-called "common envelope parameter" that describes the efficiency of deposition 
of released energy into common envelope. Trial computations have shown that the best agreement with 
observations is obtained for a cc — 2, irrespective to the evolutionary state of the donor. 

If accretor is a white dwarf, neutron star or a black hole, Eq. © is always applied. 

Nondegenerate single helium stars may be products of the evolution of close binaries. We consider three 
channels for the formation of single helium stars: the merger of two helium white dwarfs or of a helium and 
carbon-oxygen white dwarfs and the merger of a helium white dwarf and a nondegenerate helium star. In 
the first two cases the merger is a result of the loss of systemic angular momentum via gravitational waves 
radiation (GWR). Helium star and a helium white dwarf may merge either due to momentum loss via GWR 
or due to the loss of momentum in a common envelope that forms after helium exhaustion in the core of 
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Figure 1: Relation between ZAMS mass of stars M^, their mass at the end of main-sequence stage (solid 
curve), mass of helium stars produced by them (dashed curve), and mass of the final products of evolution: 
helium white dwarfs (dot-dashed curve), carbon-oxygen, oxygen-neon white dwarfs, neutron stars, and black 
holes (dotted curve) . For stars with Mi < 5 M© the upper limit of white dwarf mass attained in case C of 
mass exchange is shown. 



nondegenerate star and its expansion (for Mn e ~ (0.8 — 2.8) M Q ). If Md/M a > 2/3, the donor is disrupted 
in dynamical time scale and fo rms, most probably, a disk around its companion 

Tutu kov fc Yun gelsor Accretion in such systems is not studied as yet, but one may expect that the 

rate of accretion from the disk is defined, primarily, by the viscosity of the disk matter. The total minimum 
mass Mg" n for formation of a nondegenerate helium star by the merger of white dwarfs is, essentially, still 
a free parameter. It is clear only t hat it has to exceed the minimum mass of nondegenerate helium stars 
~ 0.3 Mq (see, for instance. Ilbenl l)l990l)). We assumed M™„ m = 0. 4 M©. The latter value is in a good 
agreement with results of computations JlbenL ll99fA lHan et all ho02h . see also discussion in Sect. The 
lifetime of merger products of helium white dwarfs or of helium stars and helium dwarfs was also calculated 
by means of Eq. though, the latter may somewhat overestimate the lifetime of merger produ cts as helium 
stars (the overestimate is < 20%, depending on the mass of the merger product, see Eq. (1) in Ilbenl (^990). 
For the merger products of helium and carbon-oxygen dwarfs, their lifetime is defined by "competition" of 
two processes: helium burning in the shell and stellar wind mass loss. We have to note that the merger of 
helium dwarfs provides in our model formation of about 90% of single helium stars. 



Based on the studies of eclipsing, spectroscopic, and visual binaries J 


Kraicheva et alJ 


1981; 


Popova et alJ 


1982; 


Kraicheva et allll985: Vereshchagin et alJ 


1988; 


Kraicheva et allll989l). we assumed that all stars are 



born in the binaries with orbital separat i ons a ranging from 6(Mi/M Q ) 1 / 3 i? to 10 6 i? Q . Birthrate function 
for binaries was taken as l|Ponova, et all 

ET 



(7) 



— =0.2d(loga) w kf(q)dq : 



(It 



where Mi is the mass of initially more massive component of the system, f(q) is the distribution of binaries 
over mass ratios of components q = M^/Mi normalized to 1. We have assumed that for close binaries 
f(q) = 1. Normalization of function JJJ assumes that one binary system with Mi > 0.8 M is born annually 
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Table 1: Birthrate and number of helium stars. 



Type of star 


Birthrate 


Total 


Subsample 




per yr 


number 


with Vn e < 16 


— c* — i 

oingle 


U.UI A 1U 


u.yu a iu 


n (\A v 1 n 4 

U.U4: A 1U 


Ivldlll - BcUUc 11 L-c LUllipdlllUll 


9i v i n -1 

U.^ A 1U 


fl 1 3 v 1 fl 7 

U.IO A _LU 


n 94 v in 3 


Xlc Wllllc-UWdll CUllipdlllUll 


nr;n v 1 n~ 2 

U.tJU A ±U 


noqy in 6 

U . <7 A _LU 


(1 IS y in 4 

U.IO A _LU 


VjW W 111 l,U 11 W CXI 1 L.U111JJ dill U 11 


n f\f{ v i n -2 

U.UU A ±U 


n i ^ v 1 n 7 

U.IO A _LU 


n fizt v 1 n 4 

U.U^fc A _LU 


( \ l\To 1X7 111 + O H TT PniTl noni Ari 
VJlNc VV 111 l,U 11 VV CXI 1 LUllipdlllUll 


nen v 1 n -3 

U.OU A ±\J 


n ^1 y in 5 

U . O -L A _LU 


(1 it; y If) 3 


Neutron star companion 


0.19 x 1(T 2 


0.31 x 10 5 


0.21 x 10 3 


Black hole companion 


0.14 x 10~ 3 


0.96 x 10 4 


0.66 x 10 2 


Binary helium star 


0.34 x 10~ 3 


0.61 x 10 4 


0.58 x 10 2 



in the Galaxy. According to 10, about 40% of all binaries are "close". Star formation rate was set constant 
for 13.5 Gyr. Evidently, this assumption does not influence the results, since the lifetime of even least 
massive helium stars does not exceed ~ lGyr. 

Our model assumes that all stars are formed in binaries. Accordingly, we did not consider the possibility 
of formation of low-mass helium stars as an outcome of evolution of single stars or of the components of 
wide binaries (see, e. g.. iD'Cruz et all l|l996ll '). The authors of the latter study assume that precursors of 
helium stars are red giants that lost almost whole hydrogen envelope immediately before helium flash or in 
the course of the flash. The mechanism of mass loss is unknown and the rate of mass loss is considered as a 
free parameter. The real rate of binarity of hot subdwarfs is difficult to determine, since for usually applied 
observational metho ds selection effects preclude detection of binar y subdwarfs with orbital periods longer 
than 200 to 300 day l|Maxted et all l200li iNa.niwotzki et all l2004al) . It is evident that a priori one cannot 
exclude the possibility of descent of low-mass helium stars from single stars or components of wide binaries, 
but, as we shall show below, it is quite possible to explain formation of sufficiently large number of single 
helium stars by the mergers of helium white dwarfs and/or nondegenerate helium stars. 

3 MAIN RESULTS 

3.1 General properties of the population of low-mass helium stars 

The main aim of our study is the modeling of the ensemble of helium stars and of their distributions over 
main parameters that, potentially, can be derived from observations - masses and orbital periods of binaries, 
as well as the analysis of the relations between parameters of close binaries with helium components. 

As it was noted in Sect. helium stars in binaries form as a result of Roche lobe overflow by stars 
that are in the hydrogen-shell burning stage. These stars may be either primary or secondary components 
of the system. In the first case the companions to the nascent helium stars are main-sequence stars, while in 
the second case they are white dwarfs, neutron stars or black holes or, on very rare occasions - also helium 
stars. 

In the Table we list total birthrate and current number of helium stars in the Galaxy, found if assump- 
tions described in Sect. are applied. The dominant constituents of the population of helium stars are 
binaries with companions - main-sequence stars and with carbon-oxygen white dwarfs. The number of the 
objects of the third main group - single stars - depends on the assumptions on the minimum mass necessary 
for He-ignition in the merger product of two helium white dwarfs: if this critical mass is reduced from 
M™ = 0.4 Mq to 0.35 Mq, the number of single helium stars increases by approximately 8%. About 0.8% 
of all helium stars have neutron star companio n s. Th e number of such binar i es may be overestimated by an 
order o f magnitude (see, e. g., iLinunov et"all l|l996h: iLhTunov et all l|l997l): IPortegies Zwart fc Yungelsonl 
(1998;i: lLommen et all 12005)), since in the present model we do not take into account the possibility that a 
nascent neutron star may get at birth a kick which may result in disruption of the system (see more detailed 
discussion in Paper II). 
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Figure 2: Distribution of helium stars over mass in the total sample (a) and in the sample limited by selection 
effects (b). Thick solid line shows helium stars with main-sequence companions, thin solid line - stars with 
white-dwarf companions, dotted line - single helium stars. 



Total model Galactic birthrate of helium stars is 0.043 yr -1 , total number of helium stars is 4 x 10 6 . 
The rate of binarity of helium stars is 76%. Since progenitors of helium stars descend from main-sequence 
stars of moderate and high mass, the semi-thickness of the Gal actic disk for them may be taken equal to 
200 p c (this agrees with empirical z-scale for helium subdwarfs l)HeberL I19861 iDownesL Il986t Ide Boer et all 
1995); then the volume of the disk component of the Galaxy is equal to 2.8 x 10 11 pc 3 , while spatial density 
of helium stars is equal to 1.4 x 10~ 5 pc ~ 3 . The latter es timate is consi derably higher than the observational 
estimates: 2 x 10~ 6 pc~ 3 (|Heherl ll98fjT . 4 x 10~ 6 pc -3 jPowneallflSfif) . However, as we shall show further, 
a considerable fraction of low-mass helium stars may be "lost" in observations due to selection effects. 

In Fig. [2] we show the distribution of helium stars with different companions over mass (for the total 
population and for the subset limited by selection effects, see below). In all groups the mass of helium stars 
strongly concentrates to the minimum value (~ 0.4 Mq). The main group of single objects consists of merger 
products of white dwarfs and white dwarfs and helium stars; very scarce single stars with mass from 2.2 
to 5.6 Mq are merger products of helium stars and helium dwarfs or former companions of binaries that 
were disrupted by supernovae explosions. In the systems with companions - helium dwarfs, the mass of 
helium stars is limited by ~ 0.9 M . This maximum corresponds to the case of completely conservative first 
mass exchange and totally nonconservative second mass exchange, when Roche lobe is filled by initially less 
massive component of the system, i. e., to the following sequence of combinations of masses of companions: 
(Mi, M 2 ) = (2.8M©,2.8M ) -> (0.4 M , 5.2 M©) -> (O.4M ,O.9M ), The upper limit of helium-star 
mass in the systems with CO- or ONe-companions is related to a similar scenario and to the upper limit 
to the mass of progenitors of white dwarfs: (Mi, M 2 ) = (11.4 M , 11.4 M ) -> (1.35 M , 21.45 M ) -> 
(1.35 Mq, 5.7 M & ). In the latter systems companion to the white dwarf explodes in the course of further 
evolution and disrupts the system. 

In Fig. 03 we show the distribution of helium stars over orbital periods. The lower limit of orbital 
periods is determined by condition for the merger of components in the common envelope. For the systems 
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Figure 3: Distribution of helium stars over orbital periods, (a) - observed subdwarfs (L. Morales- Rueda, priv. 
comm.), (b) - model sample for total population. Upper thick line shows the summary distribution, thin 
solid, upper and lower dashed lines show the input of stars with main-sequence, white-dwarf companions and 
helium secondaries into total population, (c) - distribution over periods in the sample limited by selection 
effects; the meaning of the lines is like in panel (b). 



with companions - main-sequence stars the upper limit of the periods is the maximum period that can be 
attained in the process of conservative mass exchange. For the systems with white-dwarf companions the 
upper limit of the periods is determined by reduction of orbital separations in the common envelopes. There 
are systems with P or b ~ 0.01 day in the model sample. The shortest orbital period of the known observed 
systems is P = 0.073 day (PG1017-086). Existence of model systems with shorter orbital periods may be due 
to the fact that, in the population synthesis code in the formal procedure that checks whether the merger of 
components in the common envelope is possible, the radius of a helium star with mass equal to the mass of 
the helium core of the donor is used, while a more realistic procedure has to employ much larger radius of 
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Figure 4: Distribution of companions to helium stars over mass, (a) - total model sample. Thick line - 
distribution of main-sequence companions, thin line - distribution of white-dwarf companions, (b) - model 
sample limited by selection effects; the meaning of the lines is like in panel (b). 

the inert helium core itself. It is also possible that the common envelope parameter a cc is somewhat larger 
than we assume. Another reason may be underestimate of the radii of helium stars that depend on the 
mass of hydrogen envelopes. Thus, the shortest orbital period systems shown in Fig. have, probably, to 
be classified as single helium stars. Then their number has to be lower than shown in the Figure, because 
of shorter lifetime of merger products of helium white dwarfs and helium stars. 

The distribution of the companions to helium stars over mass is shown in Fig.^t. Two discontinuities in 
the distribution for stars with main-sequence companions, occurring close to 0.3M Q and 1.0 Mq, are due to 
the loss of mass and angular momentum by main-sequence stars with 0.3 < M/Mq < 1.0 via magnetically- 
coupled stellar wind. For some systems the lifetime to contact is sufficiently shorter than the lifetime of 
helium stars with mass ~ 0.4 Mq. The upper limit of the range of mass of stars that produce magnetic 
stellar wind is still not known for certain. If this limit is raised to 1.5 M©, the number of systems with 
main-sequence companions reduces by 10%, the discontinuity shifts to 1.5 M and becomes less pronounced. 

The low-mass group of white dwarfs are helium ones, while the group of massive white dwarfs contains 
carbon-oxygen and oxygen-neon white dwarfs. Two groups of white dwarfs overlap on expense of low-mass 
(M « 0.35 — 0.40 M Q ) hybrid CO-dwarfs; precursors of the latter are helium stars of the same mass. 



3.2 Selection effects 

In the logT e ff — \ogg diagram helium remnants of stars with thin hydrogen envelopes spend their whole 
core-helium burning time in the region occupied by subdwarfs, since they have logT c ff = 4.4 — 4.5 and 
\ogg — 5.65 — 5.90. For comparison with observations, we have constructed a subset of model helium stars 
that is limited by stellar magnitude V = 1 6, since the main o bservational surveys of star s with UV-excess 
are limited by stellar magnitude B « 16 l|Oreen et all 11986) . B w 15.3 l|Downesl . Il986j) . The survey of 
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white dwarfs SPY is limited by B ss 16.5; this catalog also contains a large number of hot subdwarfs 
due to the errors of spectral classification in the input catalog. Since our model is rather primitive and 
depends on a number of not very well res tricted parameters, we neglected typical for h e lium subdwarfs 
small color excess B -V_= -0.25 ± 0.1 llHeheTl EM lAznar Cnadrado fc Jeffervl Eqq2 IStark fc Wadel 
2003: lAltmann et al.. 2004). For the systems with main-sequence companions we applied, as an effect of 
observational selection, the condition Vho < Vms- Using this rather rigid condition we get, in essence, the 
lower limit for the number of binaries where detection of companion by spectroscopic methods is possible. 
The sample limited by stellar magnitude contains only stars with mass lower than 2M , since for them 
only it appeared possi ble to construct the scale of bolometric corr ection s for V , using computations by 
llben fc Tutukovl l|l985ll . Interstellar absorption was set to l™6/kpc ijAllenl fl976f) . For all binaries we have 
assumed that for their detection the semiamplitude of radial velocity has t o be K m j n > 30 km/s; thi s 
corresponds to the lower limit of K in known binaries with helium subdwarfs JMorales-Rueda et all [2003) . 

At difference to lHan et all l|2002L 120031. we did not exclude from consideration model systems with 
spectral type G and K main-sequence companions. Han et al. sugg ested that the obje cts with Ca II 
absorption line K that were rejected as the candidates to the catalog PG 1 Green et allll98fjl . are, in reality, 



helium sub dwarfs with F- G-, and K-type main-sequence companions. However, a special analysis of rejected 
candidates {Wade et all 1200411 have shown that they really are low-metallicity sdF/sdG-subdwarfs, but not 
sdB /sdO-subdwarfs. 

Since we were interested in all helium stars, we did not take into account so called "strip effect" - 
restriction of the model sample by stars that have log g and log T e g values that are typical for core-helium 



es; in fact, some of the stars that 


Saffer et alJ 


( 


1994); 


Lisker et alJ 



The model sample of stars limited by selection effects contains 15210 objects (including systems with 
relat ivistic companions , see the Table). This number comprises ~ 0.4% of all helium stars. The catalog 
PG JGreen et allll98fj|l that is limited by B < 16 contains about 700 objects classified as helium subdwarfs 
that are distributed over the area of 10714 sq. deg. This corresponds to the total number of "observed" stars 
of about 2800. Having in mind all uncertainties of the model, rather simple account for selection effects, and 
misclassification in the PG catalog, we may consider the agreement in the numbers of stars in our model 
and in the catalog that is within factor ~ 5 as satisfactory. The most uncertain parameters of the model are 
threshold mass of the progenitors of helium stars and the minimum mass for helium ignition in the single 
merger products of white dwarfs (that constitute about a quarter of total sample) . Lowering of both these 
masses may improve agreement between the model and observations. 

Distribution of model helium stars over mass in the sample limited by selection effects is shown in 
Fig. I2t>- Reduction of the number of stars in the sample as compared to the total sample is caused, in the 
first instance, by reduction of the volume of space in which these stars might be observed (reduction factor 
~ 150 — 200). Since we exclude the systems that have main-sequence stars brighter than their companions, 
the number of helium stars with main-sequence companions further reduces by a factor ~ 20. As a result, 
in the sample limited by selection effects the rate of binarity reduces to 58%. Restriction of the sample by 
the semiamplitude of radial velocity that is necessary for discovery of binarity practically does not affect the 
number of stars in the sample. 

Observational estimates of the binary fraction of helium subdwarfs depend on the sample of stars under 
study and on the st rategy of observations, that determines the efficiency of detection of binaries. Thus, 
iMaxted et alJ J200lll estimate the binary fraction among sdB stars with periods from 0.03 to 10.0 day as 
69 ± 9%: lNapiwotzki et aTl l|2004all give (42 - 45)%, depending on as sumptions o n the r eal distribution over 
periods in the range of 0.03 to 30.0 day and masses of components: ISaffer et alJ 1120011 1 estimate the upper 
limit of the fraction of single helium stars as 35%. The reasons for the discrepancy in the observational 
estimates of binary fraction among helium subdwarfs may be differences in the samples of stars that were 
studied (they could belong to the populations wit h different age an d metallicity) and the differences in the 
observational methods. As it has been show n bvlHan et alJ feOOflh. th e fraction of binaries decreases with 
decrease of metallicity; the sample studied hv lNaniwotzki et all l|2004ajl . as admit the authors of publication, 
can contain a significant admixture of halo and thick disk stars. The model estimate of the binary fraction 
of "observed" helium stars agrees with the range of observational estimates. 

Note that reduction of common envelope parameter a cc to 1 results in increase of mergers in com- 
mon envelopes; this reduces the fraction of binaries among helium stars to 35%, in clear contradiction to 
observations. 

In the total sample of model helium stars the fraction of binaries with main-sequence components 
is 38%; selection effects reduce the fraction of such systems to 3%; these binaries may have composite 
sp ectra. The presence o f late-type companion may be evidenced by infrared excess. For instance, according 
to lStark fc Wadel {2003(1 . the fraction of systems with unresolved late-type companion among hot subdwarfs 
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may be ~ (30 — 40)%. This is close to their fraction in the total model sample. In some cases a main- 
sequence companion may be detected thanks to the reflection effect, caused by the heating by hot radiation 
of the subdwarf. But this effect may be notic e able o nly in the most close systems when companion is an 
M-type star or a brown dwarf l|Maxted at all Eq o3). Note however, that for some subdwarfs that have 
composite spectra which allow to es timate the parameters of compani ons, the masses of the latter appear to 
be confined to (0.8 - 1.2) M©range l|Aznar Cuadrado fc Jeffervl l2f)f)2j) . i. e. they are close to the maximum 
of distribution in the sample limited by selection effects (Fig. Q]d) . The rest of subdwarfs that are classified 
as single ones may be really single or have white-dwarf companions. 

While the total model population is dominated by helium stars with main-sequence companions, in the 
selection effects limited sample systems with white-dwarf companions start to dominate. The maximum of 
the distribution of helium stars over masses in total "limited" sample shifts to somewhat higher masses (into 
log(M/M Q ) = — 0.40-; — 0.35 range, see Fig- EJ. This improves the agreement with the "canonical" mass of 
subdwarfs of 0.5 M®. 

The distribution of 59 observed sdB/sdO subdwarfs with known orbital periods (Fig. Et, the data 
from L. Morales-Rueda, priv. comm.) varies within factor ~ 2 within the interval 0.1 < P or b < 10.0 day. 
In the model "observed" sample the number of stars varies within the same factor in the period range 
0.3 < P or b < 3.0 day. Relative decrease of the number of model systems to short periods is comparable to 
the observed one. With increase of P or b relative to the maximum of distribution, the number of stars in the 
model sample declines much faster than in the observed sample. There are no stars with P or b > 10.0 day in 
the latter sample. The last circumstance is possibly related to the fact that some stars which in our model 
form common envelopes in the first episode of mass exchange, in reality evolve stably to the longer periods 
(for instance, due to the loss of a part of their mass pri or to the Roche lobe overflow or due to stabilization 
of mass-exchange via momentum loss l|Han et all l200Ctl . This effect requires a test by computation of grids 
of evolutionary models. 

Accumulation of long series of observations will, most probably, result in discovery of new subdwarfs 
with long orbital periods. However, it is necessary to have in mind that the efficiency of binarity detection 
declines with increase of Por b : it exceeds 80% at P nr h < 10.0 day, decreases to 50% at P or b ~ 20.0 — 
30.0 day llMaxted et all EoOll iNaniwotzki et all l2004al) . and then rapidly declines to 0. One would expect 
that observations with high spectral resolution would result in detection of subdwarfs with relatively long 
orbital periods. But, for instance, a mong nine stars detected in SPY survey, seven have P < 1.0 day, but 
only two - P orb = 5.87 and 7.45 day l|Namwotzki et alll2004aTl . 



3.3 Relations between parameters of binaries with nondegenerate helium com- 
ponents 

In Fig. we show the number distributions of binaries with helium components and different companions 
in the selection limited model sample in the "logarithm of the orbital period log P or b - logarithm of helium 
star mass logMne" diagram. As we noted before, after taking observational selection into account, among 
"observed" helium-star binaries dominate the systems with companions - white dwarfs (~ 90%). In these 
binaries the masses of helium stars have to be confined between minimum mass of 0.35 M Q and 0.5 M®. In 
scarce systems with main-sequence components somewhat more massive helium stars may be observed, with 
masses up to (0.8 - 1.0) M®. 

Distributions of systems with different white dwarfs over periods differ: at P < 1.0 day it is more 
probable to expect the presence of a helium white dwarf than that of a carbon-oxygen or an oxygen-neon 
one. Despite we know from observations only minima of the masses of white dwarf companions to helium 
stars, let note that only in 4 out of 12 known sdB+WD systems with P or b < 1.0 day we may exclude the 
presence of helium white dwarfs (M^g 1 > 0AM®). In two known syst ems with P or b > 1.0 day whit e dwarfs 
have to belong to carbon-oxygen/oxygen-neon family (see Table 5 in llMorales-Rueda, et a,L. 2003)). 

In Fig. Elwe show, for the same model sample limited by selection effects, the distributions of systems 
with helium stars and main-sequence or white-dwarf companions in the "logarithm of the orbital period 
logP or b - logarithm of companion mass logM 2 " diagram. We plot in the Figure the masses of companions 
to helium stars known from observations, as well as the estimates of the lower limits to the masses of 
companions, if only mass function f m is known: 

= Mj sin 3 j = P olh Kf 
/m {M 1 + M 2 } 2 2vrG ' 

W e use the data on mass, semiamplitude of radial ve locity (Ki) and orbital period for 52 objects listed 
in llMorales-Rueda, et allbooaUNa.niwotzki et a,lll20f)4al) and unpublished data provided to us by L. Morales- 
Rueda. For the estimate of M™ m we assumed that Mn c =0.5 M®, sini = 1. Assumption of Mh c = 0.35 M® 
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P, day P, day 

Figure 5: Relations between orbital periods and masses of helium stars - components of binary systems that 
have secondaries of different types. The sample limited by selection effects is shown. Five grades of gray- 
scale show the regions where the number density of systems is within 1/2 of the maximum of diogpfflog m > or 
is confined to the ranges of 1/4 - 1/2, 1/8 - 1/4, 1/16 - 1/8, and - 1/16 of the maximum. Each panel is 
scaled separately. In the blank regions of the diagram systems with white dwarfs are absent. In the panels 
for binaries with white-dwarf companions, the systems that have at birth orbital periods shorter than those 
marked by the left dotted line merge in the lifetime of helium stars. Systems that have periods at birth 
shorter than the ones marked by right dotted line merge in less than Hubble time. 

only weakly influences location of stars in Fig. H3 Variation of orbital inclination from 90° to the most 
probable one ~ 60° increases M2 by a factor close to 2. In order to shift all observed objects in Fig.|J3in the 
region populated by model systems, it is necessary, as a rule, to assign to observed systems quite probable 
orbital inclination angles of 20° - 30°. 

Several binary sdB stars with M-type companions were discovered due to the eclipses or reflection effect 
in the most close systems with P < 0.25 day. As we noted above, companions with mass of 1 ± 0.2 M Q are 
suspected in some systems with composite spectra, but among systems with known orbital periods they are 
absent. 

Since the number of model stars with white-dwarf components is by a factor close to 200 higher than the 
number of systems with main-sequence companions, one may expect that the overwhelming majority of the 
new identified and discovered companions to helium subdwarfs will be carbon-oxygen or helium white dwarfs. 
As we mentioned above, the deficit of long-period model binaries may be related to the underestimate of the 
fraction of binaries with the stable first mass exchange. 

In the subsample of helium stars limited by observational selection, about 30% of stars are single. The 
mass spectrum of these stars is slightly shifted to larger periods, as compared to the spectrum of binary 
subdwarfs (Fig. EJ. These stars may be probably related to the sdB/sdO subdwarfs that have atmospheres 
enriched in He and that are slightly hotter than sdB subdwarfs with hydrogen-rich atmospheres. Location 
of these stars in the log T e g — log g diagram is more similar to the tracks of merger produc ts of helium dwarfs 
than to the tracks of single low-mass stars with helium cores and thin hydrogen envelopes l|Ahmad fc .Tef ferv. 
2003). Let note, that the o nly sample of 23 sdO-star s that was systematically explored for binarity, con- 
tains only one binary star l|Napiwotzki et all l2004ajl ! Special study of binarity detection probability by 
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Figure 6: Relation between orbital periods of binaries and masses of companions to helium stars - main- 
sequence stars and white dwarfs. The gray-scale is similar to that in Fig. In the upper panel triangles 
show known masses of main-sequence companions or their lower limits, in the lower panel triangles show 
masses or lower limits to mass of known white-dwarf companions. In both panels the crosses show lower 
limits to the mass of unknown type companions. 

iNaniwotzki et all ll20f)4al) have shown, that the binary fraction in this sample is really close to only 5%. 
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3.4 Final stages of the evolution of close binary stars with low-mass helium 
components 



As we noted above, in the closest systems with helium components the latter may fill their Roche-lobes 
before helium exhaustion in their core, because of the loss of systemic angular momentum via gravitational 
wave radiation. Respective limits of orbital periods are plotted in Fig.0for the typical mass of white dwarfs 
Af Hc = 0-2 M B Man = 0.6 M^, Mqm. = 1 . 3 M^. If conditions for stable mass loss are fulfilled, AM CVn 
systems arise l|Tutukov fc Yungelsonl Il99fil: iNelema.ns et all I2DD if) . The birthrate of suc h systems in our 
model is ~ 5 x 10 -4 yr -1 . Rather similar estimate (~ 3 x 10~ 4 yr -1 ) was obtained in iNelemans et all 
(2004) for somewhat different assumptions on star formation history and stellar evolution, but with the 
same Eq. © for description of mass exchange in the systems of main-sequence stars with comparable 
ma sses of components. Evolution of AM CVn stars an d their observed feature s are considered in detail 
by iTutukov Yungelsonl ilflQf^: INelemans et all l|200lh: INelemans et all 12004). If conditions for stable 
mass exchange are not fulfilled, the merger of helium stars with carbon-oxygen or oxygen-neon white dwarfs 
possibly results in formation of R CrB stars Jlben et all ^96) . 

After helium exhaustion in their cores, helium stars with mass < 0.8 M Q turn into carbon-oxygen or 
oxygen-neon white dwarfs of the same mass almost without any expansion. In this way close binary white 
dwarfs are formed (another scenario leads to formation of close binary white dwarfs through case C of mass 
exchange). If binary white dwarfs have short enough orbital periods, they may merge due to the angular 
momentum loss via radiation of gravitational waves in less than Hubble time. Respective limiting periods 
are shown in Fig. Currently (March 200 5) the data on orbital per iods of 19 close b inary white dwarfs is 
publis hed. Six of these systems (see, e. g., iNa.niwotzki et all l|2004hl) and Table 11 inlMorales-Rueda et all 
fcOOfil) ) have periods short enough for merger in Hubble time and one of these systems has total mass larger 
than the Chandrasekhar one. If both of com ponents of this system are carbon-oxygen white dwarfs, it may 
occur to be a precursor of type la supernova IjTutukov fc Yungelsonl Il98llllben fc Tutukovlll98l Wcbbink, 
1984). The mer ger of carbon-oxygen and helium white dw arfs, most probably, results n ot in SN la. but i n 
"helium Novae" llErgma, et aH . l200lllYoon fc La,ngeHl2004l) or in formation of R CrB stars jlhen et all ll99fi'l. 
Note, that our model predicts the rate of merger of carbon-oxygen white dwarfs with super-Chandrasekhar 
total mass equal to 0.0013 per yr , slightly lower than the prediction for th e Galactic rate of SN la equal to 
t'SNia = (0.004 ± 0.002) per yr CCapp ellaro et al l ll99flt ICannellaroL l200lll (based on the assumption that 
Milky Way has i^sNia equal to the average one for the galaxies of the same morphological type Sb-Sbc and 
blue luminosity Lb — 2.3 x 10 10 Lb,q). Theoretical a nd observational estimates of ^sNia may be brought 
to a better agreement by reduction of a ce in Eq. © llTutukov fc Yungelsonl I2002I) . O f course, it is qui te 
possible that the merger of white dwarfs is not the only mechanism for SNe la (see, e. g.. lYungelsonl ((20 05)) : 
the numerical relation between different scenarios of evolution that may potentially lead to SNe la is still 
a function of numerous parameters of the population synthesis. Note however, that our trial computations 
have shown that reduction of a ce shifts the maximum of the distribution of double helium stars to too low, 
as compared to observations, orbital periods (Fig. EJ. The situation will become more clear when sufficient 
data on close binary helium stars and white dwarfs will be accumulated. 



4 Conclusion 

We have modeled population of low-mass helium stars in the Galaxy under assumption that all of them are 
formed in close binary stars. We got estimates of their birthrate - 0.043 yr -1 , total numb er - 4x 10 6 , and bina- 
ry fra ction - 76%. These esti mates are in a goo d agreement with our analytical estimates llTutukov fc Yungelsonl 
1990) and with estimates bv lHan et all l|2003|) that are also obtained by means of population synthesis. 

Low-mass helium stars are identified with observed sdB/sdO subdwarfs. The study of observational 
selection effects shows that the dominant factor in the formation of the observed ensemble of sdB /sdO stars 
is the fact that observed samples are limited by relatively bright stellar magnitude (V < 16). The number 
of helium stars with main-sequence companions available for observations is additionally limited by the 
circumstance that majority of subdwarfs are by far brighter than their companions. All this reduces the 
binarity rate in the model observed sample to 58%; the latter value agrees with observations. According 
to our computations, in the observed population of binary helium subdwarfs overwhelming majority of 
companions to helium stars are carbon-oxygen or helium white dwarfs. 

In our model, single low-mass helium stars are, predominantly, merger products of helium white dwarfs. 
Since merge mainly low-mass dwarfs, single helium subdwarfs have masses close to those of components 
of binaries. It is possible, that a fraction of merger products may be identified with sdO subdwarfs that, 
according to existing rather scarce data, have a very low binarity rate. 

It is assumed in our model that all stars are formed in binaries, but despite this we are able to explain, 
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quite naturally, the formation of helium stars by the evolutio n in close b i naries . As well, we can explain 
the binarity rate of low-mass helium stars. On the other hand, Saffer et al. ( 20oJ), based on observations of 
radial velocities of subdwarfs, separate them into three groups. In the spectra of stars belonging to the first 
group the features of companions are absent and radial velocities of these stars are virtually constant. In 
the spectra of stars in the second group the features of companions are also absent, but variability of radial 
velocities indicates orbital periods ~ 1.0 day. In the spectra of the third group stars the features of cold 
companions are present. The percentage of stars belon ging to these three groups is 35, 45, a nd 20%. The 
following phenomenological interpretation is suggested ijSaffer et a,l.ll20oillOreen et all 1200]]^ . The stars of 
the first group are intrinsically single ob jects that lost a part of their mass not long before helium ignition 
in their cores (see scenario suggested bv iD'Cruz et al J Jl99fi^ . The stars in the second group are binaries 
that passed through the common envelope and have invisible companions - M-type stars or white dwarfs. 
In the third group there are stars that passed through the stage o f stabl e mass exchange in the red giants 
stage. If the model suggested by ISaffer et all l|200lll : lOreen et alJ ll200ll) will get additional confirmation, 
for instance, distributions of stars over orbital periods will be found, companion masses will be estimated, 
selection effects that define the ratio between different groups of stars will be established, then our model 
for the population of low-mass helium stars and the analysis of selection effects that form observed sample 
of stars would require certain revision. 

The authors acknowledge L. Morales-Rueda for providing a compilation of known orbital periods of 
subdwarfs and unpublished data. We are indebted to R. Napiwotzki, G. Nelemans, for useful remarks and 
to E. Green for providing electronic files of papers otherwise unavailable in Russia. This study was partially 
supported by Russian Foundation for Basic Research grant (project code 03-02-16254), the program of the 
Section of General Physics and Astronomy of the Russian Academy of Sciences "Extended Objects in the 
Universe", "Leading Russian Science Schools Support Program" (project code HIII-162.2003.2), and by the 
grant "Non-stationary Processes in the Universe". 
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